Nonequilibrium Phase Transition in the Kinetic Ising model: 
Dynamical symmetry breaking by randomly varying magnetic field 
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The nonequilibrium dynamic phase transition, in the two dimensional kinetic Ising model in 
presence of a randomly varying (in time but uniform in space) magnetic field, has been studied both 
by Monte Carlo simulation and by solving the mean field dynamic equation of motion for the average 
magnetisation. In both the cases, the time averaged magnetisation vanishes from a nonzero value 
depending upon the values of the width of randomly varying field and the temperature. The phase 
boundary lines are drawn in the plane formed by the width of the random field and the temperature. 
PACS number(s): 05.50.+q 



I. INTRODUCTION 

The problem of the random field Ising model (RFIM), 
having quenched random field, has been investigated 
both theoretically and experimentally || in the last 
few years because it helps to simulate many interesting 
but complicated problems. The effects of a randomly 
quenched magnetic field on the critical behaviour near 
the ferromagnetic phase transition is the special focus of 
the modern research. The recent developments in this 
field can be found in a review article ||. However, the 
dynamical aspects of the Ising system in presence of a 
randomly varying field is not yet studied thoroughly. It 
would be interesting to know if there is any dynamical 
phase transition in the presence of a randomly varying 
magnetic field. 

For completeness and continuity, it would be conve- 
nient to review briefly the previous studies on the dy- 
namic transition in the kinetic Ising model. Tome and 
Oliviera first observed and studied the dynamic tran- 
sition in the kinetic Ising model in presence of a sinu- 
soidally oscillating magnetic field. They solved the mean 
field (MF) dynamic equation of motion (for the aver- 
age magnetisation) of the kinetic Ising model in presence 
of a sinusoidally oscillating magnetic field. By defin- 
ing the order parameter as the time averaged magneti- 
sation over a full cycle of the oscillating magnetic field 
they showed that the order parameter vanishes depend- 
ing upon the value of the temperature and the amplitude 
of the oscillating field. Precisely, in the field amplitude 
and temperature plane they have drawn a phase bound- 
ary separating dynamic ordered (nonzero value of the 
order parameter) and disordered (order parameter van- 
ishes) phases. They [Q have also observed and located a 
tricritical point (TCP), [separating the nature (discontin- 
uous/continuous) of the transition] on the phase bound- 
ary line. However, such a transition, observed j?J from 
the solution of mean field dynamical equation, is not dy- 
namic in true sense. This is because, for the field am- 
plitude less than the coercive field (at temperature less 
than the transition temperature without any field), the 
response magnetisation varies periodically but asymmet- 



rically even in the zero frequency limit; the system re- 
mains locked to one well of the free energy and cannot 
go to the other one, in the absence of noise or fluctuation. 

Lo and Pelcovits [|| first attempted to study the dy- 
namic nature of this phase transition (incorporating the 
effect of fluctuation) in the kinetic Ising model by Monte 
Carlo (MC) simulation. In this case, the transition dis- 
appears in the zero frequency limit; due to the presence 
of fluctuations, the magnetisation flips to the direction 
of the magnetic field and the dynamic order parame- 
ter (time averaged magnetisation) vanishes. However, 
they (|] have not reported any precise phase bound- 
ary. Acharyya and Chakrabarti ^ studied the nonequi- 
librium dynamic phase transition in the kinetic Ising 
model in presence of oscillating magnetic field by exten- 
sive MC simulation. They |J have successfully drawn 
the phase boundary for the dynamic transition and ob- 
served/located a tricritical point on it. It was also no- 
ticed by them Q that this dynamic phase transition is 
associated with the breaking of the symmetry of the dy- 
namic hysteresis (m — h) loop. In the dynamically dis- 
ordered (value of order parameter vanishes) phase the 
corresponding hysteresis loop is symmetric, and loses its 
symmetry in the ordered phase (givingnonzero value of 
dynamic order parameter). They have [|9| also studied the 
temperature variation of the AC susceptibility compo- 
nents near the dynamic transition point. They observed 
that the imaginary (real) part of the AC susceptibility 
gives a peak (dip) near the dynamic transition point 
(where the dynamic order parameter vanishes). They 
H have the following conclusions: (i) this is a distinct 
signal of a phase transition and (ii) this is an indication 
of the thermodynamic nature of the phase transition. 

Recently, the relaxation behaviour, of the dynamic or- 
der parameter, near the transition point, has been stud- 
ied both by MC simulation and solving meanfield dy- 
namic equation. It has been observed that the relaxation 
is Debye type and the relaxation time diverges near the 
transition point, showing a critical slowing down. The 
'specific heat' and the 'susceptibility' also diverge |Tl| ] 
near the transition point in a manner similar to that of 
fluctuations of the dynamic order parameter and energy 
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respectively. 

The statistical distribution of the dynamic order pa- 
rameter, near the dynamic transition point, has been 
studied by Sides et al E^] . They observed that the distri- 
bution widens up (to a double hump from a single hump 
type) as one crosses the dynamic transition phase bound- 
ary. Since, the fluctuation increases as the width of a 
distribution increases, this observation is consistent with 
the independent observation |lf[] of critical fluctuations of 
dynamic order parameter. They [[12[ have also observed 
that the fluctuation of the hysteresis loop area becomes 
considerably large near the dynamic transition point. 

In this paper, the dynamic phase transition has been 
studied, in the two dimensional kinetic Ising model in 
the presence of randomly varying (in time but uniform 
in space) magnetic field, both by MC simulation and by 
solving the meanfield dynamical equation. This paper 
is organised as follows: in section II the model and the 
MC simulation scheme with the results are given, in sec- 
tion III the meanfield dynamical equation and its solution 
with the numerical results are given and the paper ends 
with a summary of the work in section IV. 



II. MONTE CARLO STUDY 

A. The Model and the simulation scheme 

The Hamiltonian, of an Ising model (with ferromag- 
netic nearest neighbour interaction) in presence of a time 
varying magnetic field, can be written as 



H 



<ij> 



J 11 & 4 ' 



(2.1) 



Here, s\ (= ±1) is Ising spin variable, Jjj is the interac- 
tion strength and h(t) is the randomly varying (in time 
but uniform in space) magnetic field. The time variation 
of h(t) can be expressed as 



h{t) 



h r(t) 




for 
elsewhere 



t Q < t < t + T 



(2.2) 



where r(t) is a random variable distributed uniformly 
between -1/2 and +1/2. The field h(t) varies randomly 
from — ho/ 2 to ho/ 2 and 



t +r 



h(t)dt = 0. 



(2.3) 



The system is in contact with an isothermal heat bath 
at temperature T. For simplicity the values of all Jy are 
taken equal to unity. The periodic boundary condition is 
used here. 

A square lattice of linear size L(= 100) has been con- 
sidered. Initially all spins are taken directed upward and 
h(t) = 0. At any finite temperature T, the dynamics of 



this system has been studied here by Monte Carlo sim- 
ulation using Metropolis single spin- flip dynamics. The 
transition rate (sf — > — sf ) is specified as 



W(sf 



-sf ) = Min [1, exp(-AH/K B T)} (2.4) 



where AH is the change in energy due to spin flip and Kb 
is the Boltzmann constant which has been taken unity 
here for simplicity. Each lattice site is updated here se- 
quentially and one such full scan over the lattice is defined 
as the time unit (Monte Carlo step per spin or MCSS) 
here. The magnitude of t he fi eld h(t) changes after every 
MCSS obeying equation f2~2[ The instantaneous mag- 
netisation (per site), m(t) — (1/£ 2 )X^ S ? nas been cal- 
culated. After bringing the system in steady state (m(t) 
gets stablised with some fluctuations), the switch of the 
randomly varying magnetic field h(t) has been turned on 
(at time to MCSS) and calculated the instanteneous mag- 
netisation. The to has been taken equal to 2 x 10 6 and 
even more (3.15 x 10 6 ) near the static ferro-para transi- 
tion temperature (~ 2.269..). By inspecting the data and 
the time variation, it is observed that m(t) gets stabilised, 
for this choice of the value of to- The time averaged (over 
the active period of the magnetic field) magnetisation 
Q = 7 Jt^ +T rn(t)dt has been calculated over sufficiently 
large time r (1.75 x 10 6 ). By changing the values of r 
the stabilisation of Q has been checked quite carefully. 
It is observed that, Q does not change much (apart from 
the small fluctuations) and this choice of the value of r 
is considered to be good enough to believe that the value 
of Q becomes stable. Here, Q plays the role of dynamic 
order parameter. It may be noted that, one measures 
the value of order parameter i.e., sponteneous magneti- 
sation, in the same way in the case of normal ferro-para 
transition. But, there is a remarkable difference. In that 
case, system reaches a steady equilibrium state, however 
in this case the state lies in a nonequilibrium (time depen- 
dent) state. This dynamic order parameter Q is observed 
to be a function of width (h ) of the randomly varying 
magnetic field and the temperature (T) of the system, 
i.e., Q = Q(ho,T). Each value of Q has been calculated 
by averaging over at least 25 different random samples. 



B. Results 

Taking all spins arc up (sf = +1) as the initial condi- 
tion, the simulation was performed using the above form 
of the time varying field h(t). It has been observed nu- 
merically that, for a fixed values of ho, if T increases, Q 
decreases continuously and ultimately vanishes at a fixed 
value of T. Similarly, for a fixed temperature, if ho in- 
creases, the value of Q decreases and finally vanishes at a 
particular value of ho- Figure 1 shows the time variation 
of magnetisation m(t) at a particular temperature T and 
for two different values of field width ho. For a small 
value of ho, the system remains in a dynamically sym- 
metry broken phase (Fig. la) where the magnetisation 
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oscillates asymmetrically about zero line. As a result, the 
dynamical order parameter Q (the time averaged mag- 
netisation) is nonzero. As the held increases, the system 
gets sufficient energy to flip dynamically in such a way 
that the magnetisation oscillates symmetrically (Fig. lb) 
about zero line and as a result Q vanishes. 

It is possible to let Q vanish, either by increasing the 
temperature T for a fixed held width ho or vice versa. 
It has been observed that for any fixed value of ho the 
transition is continuous. Two such transitions (for two 
different values of ho) is depicted in Fig. 2. So, in the 
plane formed by the temperature (T) and the held width 
{ho), one can think of a boundary line, below which Q 
is nonzero and above which it vanishes. Figure 3 dis- 
plays such a phase boundary in the ho—T plane obtained 
numerically by Monte Carlo simulation. The transition 
observed here is continuous irrespective of the values of 
ho and T, i.e., unlike the earlier cases (^^], no tricritical 
point has been observed. It may be noted here, that the 
transition temperature, for held width going to zero, re- 
duces to the equilibrium (zero held) ferro-para transition 
point. In hg. 3, the overestimated (from the Onsager 
value T c = 2.269..) value of the transition temperature 
in ho — > limit, may be due to the small size of the sys- 
tem. It should be mentioned here, that all the results 
are obtained by using sequential updating scheme. Al- 
though, there is some study fl3|| regarding the updating 
scheme to simulate the dynamic processes, in the present 
study no significant deviation was found for these two 
different (sequential/random) updating techniques. Few 
data points of the phase boundary (in Fig. 3 marked 
by bullets) are obtained by using the random updating 
scheme and no significant deviation was observed. 



III. MEANFIELD STUDY 

A. Meanfleld dynamical equation and solution 

The meanheld dynamical equation of motion for the 
average magnetisation (m) is 



dm 
~~dt 



-m + tanh 



m{t) + h{t) 
T 



(3.1) 



where h{t) is the randomly varying external magnetic 
field satisfying the following condition 



1 



to+r 



h(t)dt = 0, 



(3.2) 



to 



with the same distr ibution P{h) discussed earlier. 

The equation 3.1 has been solved emplyoing the fourth 
order Runge-Kutta method subjected to the above con- 
dition. The initial magnetisation is set equal to unity, 
which serves as a boundary condition. First, the mag- 
netisation m{t) has been calculated for h{t) = and 



brought the system into an equilibrium state and af- 
ter that the switch of the randomly varying magnetic 
field has been turned on. The dynamic order parameter 
Q = i Jq m{t')d t' ha s been calculated, from the solution 



m{t) of equation 
random samples 



3.1 



The Q is averaged over 25 different 



B. Results 

Observations similar to the MC case are made in this 
case. For quite small values of ho and T systems remains 
in a dynamically asymmetric phase {Q ^ 0) and gets into 
a dynamically symmetric {Q = 0) phase for higher values 
of ho and T. It is observed that one can force Q to van- 
ish continuously by tuning ho and T. A phase boundary 
line for the dynamic transition is shown in Fig. 4. Here 
also, the limiting {ho — > 0) trantistion temperature re- 
duces to the MF equilibrium (zero field) transition point 
(j>mf = iy it may be noted here that the coordination 
number z {= 4 in two dimension) has been absorbed in 
the interaction strength J in calculating the T^ IF {= 1) 
and the temperature (shown here in Fig. 4) is measured 
in the unit of Jz. 



IV. SUMMARY 

The nonequlibrium dynamic phase transition, in the 
two dimensional kinetic Ising model in presence of a ran- 
domly varying (in time but uniform over the space) mag- 
netic held, is studied both by Monte Carlo simulation 
and by solving the meanfield dynamical equation of mo- 
tion. In both the cases, it is observed that the system 
remains in a dynamically symmetric phase {Q = 0) for 
large values of the width of the randomly varying mag- 
netic field and the temperature. By reducing the value of 
field width and temperature one can bring the system in a 
dynamically symmetry broken phase (Q ^ 0). The time 
averaged magnetisation, i.e., the dynamic order param- 
eter vanishes continuously depending upon the value of 
width of the randomly varying field and the temperature. 
The nature of the transition observed here is always con- 
tinuous and the phase boundaries are drawn. It may be 
mentioned here that the dynamic responses, of Glauber 
kinetic Ising model, are studied || for a quasiperiodic 
time variation of the magnetic field. 

There is some experimental evidence of dynamic tran- 
sition. Recently, Jiang et al jl4| observed the indica- 
tions of dynamic transition, associated with the dynam- 
ical symmetry breaking, in the ultrathin Co/Cu(001) 
sample put in a sinusoidally oscillating magnetic field by 
magneto-optic Kerr effect. For small values of the ampli- 
tude of the oscillating field the m — h loop lies asymmet- 
rically (dynamic ordered phase) in the upper half plane 
and becomes symmetric (dynamic disordered phase) for 
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higher values of the field amplitude. Very recently the 
dynamical symmetry breaking has been observed experi- 
mentally in highly anisotropic (Ising like) ultrathin ferro- 
magnetic Fe/W(110) films and nicely depicted in Fig. 1 
of Ref . . However, the detailed quantitative study to 
draw the phase boundary is not yet done experimentally. 

Very recently, the stationary properties of the Ising 
ferromagnet in presence of a randomly varying (having 
bimodal distribution) magnetic field are studied Jl6| in 
the meanfield approximation. The transition observed 
from the distribution of stationary magnetization is dis- 
continuous. However this present study deals with the 
dynamical properties (dynamical symmetry breaking) of 
the system and the noncquilibrium transitions in the 
Ising ferromagnet in presence of randomly varying (uni- 
formly distributed) magnetic field. In the former case the 
transition was observed from the distribution of station- 
ary magnetisation. However, in the present study, this 
was observed from the temperature variation of dynamic 
order parameter associated with a dynamical symmetry 
breaking and the transition observed is continuous. An 
extensive numerical effort is required to characterise and 
to know the details of this dynamical phase transition. 
It would also be important to see whether the different 
kind of distribution of the randomly varying field would 
give different results. 
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Figure captions 

Fig. 1. Monte Carlo results of the time (t) variations 
of randomly varying field h(t) and the response magneti- 
sation m(t) for T = 1.7 and (a) ho = 1.0 and (b) ho 
TO. The symmetry (about the zero line) breaking is clear 
from the figure. 

Fig. 2. Monte Carlo results of the temperature (T) 
variations of the dynamic order parameter Q for two dif- 
ferent values of ho. (•) for ho = 2.4 and (O) for ho = 
0.8. 
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Fig. 3. The phase boundary for the dynamic transition 
obtained from MC simulations. The data points repre- 
sented by (•) are obtained by using random updating 
and those represented by (O) are obtained by sequential 
updating. 

Fig. 4. The phase boundary for the dynamic transition 
obtained from the solution of the MF dynamic equation 
(equation 3.1). 
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